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DDttDOM 


A MODIFIED REYNOLDS ANALOGY FOR THE COMERESSIELE 
TURBULENT BOUNDARY LAYER ON A FLAT PLATE 
By Morris W. Rubesin 


SUMMARY 


A modified Reynolds analogy is developed for the compressible 
turbulent boundary layer on a flat plate. When mixing -length theories 
are used to evaluate terms of the final expressions, it is found for 
air that the ratio of Stanton number to half the local skin-friction 
coefficient is greater than unity. At Mach number equals zero, this 
ratio is of the order of I .18 to 1.21 for Reynolds numbers based on 
momentum thickness of 10^ to 10^. Up to a Mach number of 5 snd under 
extreme conditions of surface temperature, it is found that the ratio 
of Stanton number to half the skin-friction coefficient differs from its 
values for the inccmpressible case (M=0) by amounts so smRn as to be of 
the magnitude of the uncertainties in the theory. 


INTRODUCTION 


There are several theories in the literature of aerodynamics which 
are concerned with the subjects of skin friction and heat transfer in 
the compressible tvirbulent boundary layer on a flat plate (refs. 1 
through 8 ). Each of these theories, however, is restricted throu^ the 
assumptions that Preindtl number (Pr) is unity in the 1 B.im‘ uar sublayer 
and that there is an equivalence in the mechanisms of the transport of 
heat and momentum in the turbulent region of the bo-undary layer. ^ The 
latter can be considered equal to the assumption that the turbulent 
Prandtl number (a) is unity. From these assiimptions and from a definition 
of the heat -transfer coefficient based on the temperature difference 
between that of the siirface and the stagnation temperatirre of the free 
stream, it is foimd that there is an exact equivalence bet^reen the local 

heat-transfer coefficient (written in dimensionless fashion as the 

iThe recent mixing-length theory of Lin and Shen (refs. 6 , 7^ and S) 
considers the differences in the turbulent exchange mechanisms of 
momentum and energy. The theory, however, is incomplete in that it 
requires empirical knowledge of three coefficients which have not yet 
been determined.. 
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Stanton number) and the local skin-friction coefficient. 03iis corre- 
sponds to tbe -weU-known Reynolds analogy of the low-speed case. 

Rrom investigations of beat transfer at subsonic speeds ^ it is 
known, for the case of air, that the Reynolds analogy underestimates the 
value of the Stanton number by approximately 20 percent. The results 
from these investigations, including the enpirical work of Colburn 
(ref. 9 ), the results of this analysis are shown in table I 
(refs. 10 t.h-rn ugh l4). The analyses in the table are classified 
according to (l) which of the usual three boundary-layer subdivisions 
are used, (2) what assTJnptions are employed in each, (3) what value 
of turbulent Prandtl number is enployed, and (4) whether frictional 
dissipation is included. In general, the following is observed; 

Including the buffer layer iptroduces to the final expressions an 
additional tenii which depends on the ass^^med velocity distribution in 
that region; including dissipation leads to the evaluation of a tenper- 
ature recovery factor; and allowing for variation of the turbulent 
Prandtl number a results in terms of Pr/a where Pr appeared for 
the cases where a = 1. It should be emphasized that only Shlrokow, 

Smith and Earrop, and the present analysis consider the effect of 
conpressibility. The other analyses neglect compressibility and employ 
low-speed values for the velocity ratios at the outer edges of the 
la m i n a r sublayer and buffer layer, respectively, and use empirically 
determined inconpressible velocity distributions in evaluating the 
effect of the buffer layer. The latter procedure introduces the constants 
which appear in the end expressions developed by von K&m^ and Seban. 

It should be noted that Shlrokow also restricts his end results to the 
incompressible case by employing a low-speed veilue for the velocity 
ratio at the outer edge of the sublayer. 

The question arises as to how the relative thickening of the sub- 
layer and buffer layer under the conditions of compressibility will 
Influence the results of these analyses. Will compressibility caxxse 
the large corrections to the Reynolds analogy to increase or decrease? 
Smith and Harrop (ref. l4) have examined this problem from the compressi- 
bility viewpoint. The identical assumptions used by von TTa-rmnn were 
employed except that all temperatures were replaced by total temperatures 
to account for frictional dissipation. This procedure is justified 
whenever the lami nar or tvcrbulent Prandtl nimiber is unity; as Smith and 
Harrop assumed for the turbulent portion of the boundary layer. Since 
the analysis considered the laminar Prandtl number other than unity in 
the sublayer and buffer layer, seme error was introduced. The mn-tn 
effect of this error was that the recovery ten^rerature of the surface 
remained at the stagnation temperature whatever the value of Pr. In 
general, when it was assumed that the surface temperature largely 
governs the extent of the sublayer and buffer layer, it was found that 
compressibility has a very sm all effect on the relationship between heat 
transfer and skin friction. 
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In view of the failure of the Smith and Haxrop theory to show any 
■variation of tanperature recovery factor -with Prandtl number and of the 
omission of the effect of the turbulent Prandtl number, it is believed 
desirable to re-examine the problem considering -the effects ingposed by 
compressibility. It is the purpose of this report, therefore, to study 
■the relationship between heat transfer and sMn friction in the com- 
pressible turb'ulent boundary layer on a flat plate ■with en^hasis on 
the following; 

lo The effect of "the relati"va thickening of the sublayer 
of the boundary layer 

2 . The effect of frictional dissipation 

3. The effect of differences in the turbulent exchange 

mechanism of momen'tum and heat (turbulent Prandtl 
ntimbers other than unity) 

Because of the present lack of unders^tanding of the mechanism of 
turbulence in shear flow and of -the structure of compressible turbulent 
boundary layers, ■this theory, like all ■fehe o^ther theories of the 
turbulent boundary layer, depends largely on arbitrary, but plausible, 
assumptions. It cannot, ■therefore, be considered as absolute. 


SYMBOLS 


Cf 


E 

h 

H 

k 

M 


local skin-friction coefficient. 


(l/2)p^u„2 


specific heat at cons^bant pressure 
to^bal energy, CpT + u^/2 

loceil heat-transfer coefficient, rate of heat ■transfer 
per unit area per degree, q7(T^Tav) 

thermal energy flux (defined by eg.. (8)) 

thermal conducti^vl^ty 

Mach number 


p pressure 

Pr Prandtl number conrposed of properties based on molecular 
transport, UCp/ij; 

q. rate of heat transfer per unit area 
Rq Reynolds number based on momen^tum thickness 
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r local recovery factor, 

St Stanton number, b/p^u^Cp^j^ 

T temperature 

u velocity in the x direction 
V velocity in the y direction 

x,y coordinate system with y = 0 plane as surface of the 
plate 

a Prandtl n umb er coirposed of properties based on turbulent 
transport, eCp/u 

7 ratio of specific heats ( 1.4 for air) 
e eddy viscosity (defined by eq.. (9) ) 

0 momentum thickness 

K eddy thermal conductivity (defined by eq. (lO)) 
p viscosity 

p density 

T shear stress (defined by eq. ( 7 )) 

Siiperscripts 

— mean value 

’ fluctuating value 

Subscripts 

1 condition at outer edge of sublayer 

2 condition at outer edge of buffer layer 

S stagnation conditions 

condition at outer edge of boundary layer 
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w condition of sxirface 

aw condition at surface for zero heat transfer 

System of units used is arbitrary .but must be self consistent . 


AMLYSIS 


The basic compressible turbulent botmdary-layer equations are 
given in reference 5 can be written for steady flow past a flat 
plate as 




(puv) = (u ^ - u p’v* 
oy dy 


p u'v’ ) 


( 1 ) 


(p{lE) + 1- (pvE) = A (- g _ p 




ay Sy 


E'v' - E p’ v' )+ 




( 2 ) 


^ (pu) + ^ (pv + p'v») = 0 


(3) 


At present, equiations (l) to (3) cannot be solved rigorously. They 
can only be used as guides toward choosing the Important variables and 
for indicating the form of simplified, thou^ arbitrary, relationships 
between these variables. 


For choosing the variables to be considered, equations (l) to ( 3 ) 
are rewritten using the relation 


E = cpT + |- u2 


i^) 


as 


— aii 

pu 


^ + (pv + p>v') ^ ^ (fl ^ - p u'V') 


ay ay ay 


(5) 


^ (k ^ + ujl p CpV’T' - up u'v') 


ay ay 


ay 


(6) 
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In equations ( 5 ) and ( 6 ), the fluctuating density terms have been 
placed on the left side of the equations. By virtue of the continuity 
equation ( 3 ), It Is plausible to consider the fluctmtlng density terms 
as virtviELl masses which contribute to those terms -vdiich represent the 
change of momentum or of total energy of a unit volume of fluid at a 
point. As a resxalt of this, the fluctuating density terms do not 
appear in the terms which are considered to represent the transport of 
momentum and of thermal energy on the right of the equations. 

The terms within the operators on the ri^t side of equations ( 5 ) 
and (6), called the shear stress and energy flux, are defined as 

T = A ^ - p u’v' = (u + e) ^ (7) 


H =£ 


— ^ 


- p Cp v’T’ + up ^ 


- up u*v* 


= (k+ k) ^ + u(ii + e) 


where the eddy viscosity e is defined as 


p u*v* 


( 8 ) 


(9) 


and the eddy thermal conductiviuy K is defined as 


p Cp v'T' 

. " - - (1°) 

The bars have been dropped in the final terms of equations (7) and (8) 
because all the terms represent mean values. 

Although equations ( 5 ) and ( 6 ) cannot be solved rigor oiisly, they 
do yield sufficient information to act as a guide in treating the 
problem approximately. For the case where Pr = 1 and a = 1, it can 
be shown that the dependent variables in equations ( 5 ) and (6) are 
linearly related to each other, that is, 

CpT + ^ = Au + B 

idiere A and B are constants. This - is equivalent to expressing H/t 
equal to a constant through the boundary layer. For the case where 
both Pr pirirl a differ from but are close to unity, it is plausible to 
assume, as a first app3roximation, that H/t still rema ins constant 
through the boundary layer, thou^ the constant may be dependent 
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on Pr and a. It should be noted that this assumption, that h/t = con- 
stant, is the basis for the present theory. 

The ratio of the terms of equation (8) to those of equation (7) 
is 


1 = k + B . ^ 
|I .+ € 8u 


+ u 

X 


( 11 ) 


The relation between temperature and velocity in the boundary layer is 
obtained on integration of equation (ll) 


T - 




+ € 


(| - u) du 


( 12 ) 


At the surface of the plate, the molecular treinsport terms pre- 
dominate and equations (7) and (8) become 


Therefore, 



( 13 ) 


Equation (12) becomes 

Equation (l4) can be rewitten as 



pU 



1 

1 

II 

u + e , 

'*’w 

Jo k + K 

du 


(14) 


( 15 ) 


For the case of zero convective heat transfer at the surface, the 
numerator of equation (l5) nrust be zero, requiring the surface temper- 
ature to be given by 


(1 + 




-P- + g 

k + K 


Tw — Taw ~ T 


u du) 


(16) 
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or 


■-aw 


= 





1^ + e 
k + K 



(17) 


The quantity Tav is called the recovery tengjerature . When the 
heat-transfer coefficient is defined as 


q.w 


(18) 


the ratio of Stanton number to skin -friction coefficient becomes 

h 


St 


_ Pg^ca ^P«P _ 


Cf 

2 


2 k + K 


From equation (16) and the definition of recovery factor 

Taw “ T 


r = 


TS - T„ 

the e3q)resslon for recovery factor is 


(19) 


( 20 ) 


r 



Cpoo (1^ +e) 

k + K 



s 


( 21 ) 


Equations (l 9 ) and ( 21 ) constitute the general expressions for the 
ratio St/(cf/ 2 ) and the recovery factor r. Integration of these 
equations is facilitated by the assumption that the boundary layer can 
be divided into three parts: the sublayer, where only the molecular 

trainsport terms appear (integrand = Er); the turbulent portion, where 
only the eddy transport terms appear (integrand = <x); and the buffer 
lajer, ■\daere both kinds of transport occur, the integrand depen ding 
on the proportion of each type of transport mechanism. To establish 
the proportion of turbulent to molecular exchange in the buffer layer, 
von Karm^, Reichardt, and Seban used velocity distributions to estab- 
lish e and then confuted tc by making an assumption of the value 
of a. (See table I.) Smith and Harrop carried this over to the com- 
pressible case by assuming that the velocity distribution in the buffer 
layer i-Tas Identical to that used by von Karm^ when the fluid properties 
are evaluated at the temperature and when the velocity term is 

arbitrarily replaced by a function of velocity obtained from their 
expressions for the fully turbulent-flow region. Because of the 
uncertainty of the latter procedure, the equally uncertain but simpler 
approach of entirely neglecting the buffer region and some'VThat thicken- 
ing the sublayer will be used in this analysis. It can be deduced from 
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table I by comparing the Prandtl-0?aylor theory "with von Ka rniRn » r theory, 
for the case of O .5 < Er < 2, that a very small error is introduced by 
omitting the buffer region and replacing it by a thicker sublayer. 

When the integrand in equations ( 19 ) and (21) is set equal to Pr 
in the sublayer, and a in the turbulent portion, the equations became 

St ^ 1 

£f %[i- 

and 

r = a [1 - (1 - J 

*-*00 

It is observed from table I that for a = 1, equation (22) is 
identical with those derived by Taylor, Erandtl, ani^ Shirokow. For 
a = 1 , the recovery-factor expression is identical with Shirokow' s- 

When a is eliminated between equations (22) anfl (23) there 
results 


( 22 ) 


(23) 


Ui 

St ^ ^ ^ 

^ r + Er 

o 


(24) 


DISCUSSION 


The results of this analysis have been shown to depend on a knowl- 
edge of the value of the velocity ratio at the outer edge of the sub- 
layer and of the value of the turbulent Erandtl number a. Neither of 
these quantities can be expressed with certainty; however, it will be 
shown that siofficient information is available to allow numerical 
evaluation of the ratio of Stanton number to the local skin-friction 
coefficient . 

There are alternative procedures for the evaliiation of u^/u^ . 

The dost widely used method is to assume that low-speed expressions 
prevail \mder the conditions of campressibillty with the exception that 
the fluid properties are evaluated at the surface tenperatinre . Accord- 
ingly, Frankl and Voishel (ref. 1) express this as 



(25) 
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The use of this relationship in Frank! and Voishel’s theory for 
skin friction (also in^jlied in Van Driest 's theory) is 1010 ™ to yield 
ski n -friction coefficients ■which agree irell "with a-vailable experimental 
data at Mach numbers up to 2.5 (ref. 15) • Equation ( 25 ) 'will be used 
in the nmerical ■e'valuations "which follow; however, effects of the 
uncertainties introduced •will be noted. 

In order to use equation ( 25 ), it is necessary to be able to 
e'valuate the local skin -friction coefficient cf . At present there are 
insufficient experimental data to establish cp under general conditions 
of surface temperature and Mach number. It is necessary to rely on 
theory. For this analysis the results of the appendix, in which the 
well-known Van Driest analysis is repeated for Er / 1 and a 5 ^ 1, are 
■used in e~yaluation of cp. 

Fresent knowledge of the turbulent Frandtx numoer a is also quite 
limited. The term a, -unlike Pr, is believed to depend on the type of 
flow as -well as on the fluid because differences have been found in low- 
speed tests of jets and of pipes or channels. In air jets it has been 
found (ref. 1(5) that a -varies from O .70 to 0.77* These values of a 
are over-all -values determined from the spreading of ’the jets. In a 
pipe it has been found (ref. 17 ) that the local a varies from 0.82 
to 1 . 06 , depending on the radial position and on the Reynolds number. 
Experiments in rectangular channels (ref. I 8 ) also indicate a variation 
of locEil a ranging from O .60 to 0.95 due to position and Reynolds 
number. There are no data a-vailable, however, for -the magnitude of a 
in boundary layers on a flat plate. 

In -view of -the . arbitrary character of -the assumption concer ning the 
velocity ratio at the outer edge of -the sublayer and of the uncertainty 
of the value of gt, can the res-ults of -this theory yield information 
of any accuracy? It is believed that equation (24) allows a fairly 
accvirate estimation of ‘the ratio St/(cp/2). The reasoning -which leads 
to this concltislon is based on -fche considerations of "the following 
paragraph. 

When a constant value of a, having the same order of magnitude 
as is measured subsonically in pipes, and the value of u^^/ujo based on 
equation ( 25 ) are in-troduced into equation ( 23 ), it is found that the 
recovery factor is louvered, from its incompressible -value by as much as 
10, percent at a Ifech number of 4. This is contrary to experimental 
experience where the recovery factor r remains essentially constant 
■with Mach number in the range 0 < M < 3*8 (refs. I 9 , 20, 21, and 22). 

To conform -with experiment, a or ux/u„ or both would have to be altered 
as a function of Mach number to maintain the recovery factor constant. 

To avoid this speculation, the direct influence of a is eliminated 
through -the use of equation (24) where the experimentally verified 
expression r = Rr^/3 is -used. The primary source of error in 
equation (24), -then, is -the -uncertainties in the value of 
Because of "the nature of ecguation (24) "tdien Pr is near 0.7, i'ts 



NACA TN 2917 


11 


approximate value for gases, uncertainties in ux/uoo are not too 
serious. On the average, in the range 0 < M < 5 a 20-percent error 
in Ux/uja represents a 1-percent eiror in the ratio St/(cf/2). 

Values of St/(cf/2) determ in ed from equation (2^) are shown in 
figures 1 and 2 for Pr = O. 72 . The two limiting surface temperatures 
likely to be encountered are shown, respectively, in the two figures. 
These are the recovery temperature and the free-stream teurperature . For 
the case of M = 0 both surface temperatures are equal, and the results 
of figures ,1 and 2 indicate that the ratio of the Stanton number to 
half the local skin-friction coefficient is from about 17*5 percent 
to 20.5 percent higher than the Reynolds analogy. Colbiim's (ref. 9) 
empirical correction to Reynolds analogy is Pr-2/3 and corresponds to 
a 25 -percent correction. The enq)irlcal correction, however, is based 
on data below a length Reynolds number of one million where the correc- 
tion of this theory would be about 21 percent. This correspondence of 
theory a n d experiment is well within the scatter of the experimental 
data. For the case of a surface at very near the recovery temperature, 
it is noted from figure 1 that the Influence of Mach number is small, 
increasing the ratio of the Stanton number to half the local skin- 
friction coefficient at M = 5 ^y only about 3 percent over the incom- 
pressible case. For the cooled case the effect of Mach number is n gn^n 
small, lowering the ratio St/(cf/2) by about 1 percent at M = 5 . 

These variations can be considered to be the same order of magnitude as 
the uncertainties in the theory for Rrandtl ninnber near unity. 

A comparison of the results of the present theory with the results 
of Smith and Harrop for Rg = 10® is also shown in figure 1. It is 
observed that the Smith and Earrop theory yields results which are 
markedly lower than the results of the present theory. By varying the 
reference velocity for the buffer layer (uq in table I) between its 
limits, it was found that the Smith and iferrop values varied by less than 
3 percent. The main difference in the results of the two theories, 
therefore, is due to the consideration of recovery factor in the present 
theory. A similar comparison to that shown in figure 1 cannot be made 
in figure 2. For the condition Tw = T^o, it is found that the series 
in the denominator ia Smith and Earrop ’s end equation (table I) does 
not converge rapidly enough to give a valid answer at Mach numbers of 
the order of luaity or higher. 


CORCLiroiRG REMARKS 


From a simple analysis of the effect of compressibility on the 
relationship between heat transfer and skin friction for air, a modified 
Reynolds analogy, it is found that the ratio of Stanton number to 
the local skin-friction coefficient at M = 0 is from 17 *5 percent 
to 20.5 percent hi^er than given by Reynolds analogy in the range of 
momentimi-thickness Reynolds number from 10® to 10®. 
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It is also found that up to a Mach nujnber of 5 under extreme 
conditions of surface ten^jeratures , the effects of compressibility on 
this ratio are so small as to be of the same order of magnitude as the 
possible uncertainties of the theory. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Dec. 8, 1952- 
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APEEMDIX 

SKIN-FRICTION ANALYSIS 


In this appendix, the Van Driest analysis will be repeated with the 
modifications that Pr 7 ^ 1 and, a 1 . 

The basic equation for determining the velocity distribution in the 
turbulent portion is 

T = p E2y2 (Al) 

where the quantj.ties represented are temporal mean values. 

By algebraic manipulation, equation (Al) is transformed to 



(A2) 


Equation (A2) corresponds to equation (48) of reference 2. Letting 
iy = u/u^ equation (A2) becomes 


d\!C _ 1 /tv dy 


(A3) 


Since the static pressure across the boundary layer is considered 
constant 


Pw _ T 
P " 


m 


On integration of equation (12) of the text, the ri^t member of 
equation (A4) can be written as 


|- = b2 (1 + - A2 u2) 


where 
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= 1 - nip (1 - fe) + n^p2(l - fl) 

1 + 2^ a [ 1 - (1 - ^)^^]- ^ 

^ [1 - (1 - ^)p] 


n2 


7-1 


a M, 


Tw 


T 

■^00 




(a6) 


Men a = Pr = Ij eqmtions (A5) and (a 6) degenerate to those used by 
Van Driest. Substitution of equation (A5) into equation (A3) resiolts 
in 


du _ b 

y 1 + bOl - A 25 T 2 KU 


dy 

V y 


(AT) 


Except for the term b, equation (AT) is identical with equation (52) 
of Van Driest. Following Van Driest, equation (AT) is integrated to 
yield 



2A^ u - B 
(b 2 + li-A^) V2 


+ lain-i 


B 

(b^ + 4a^)^/^ 



where F is a constant. 

Men equations (a 4), (A5), and (a 8) are substituted into 


there results 




- u) dy 


^Pw d 
S dx 


[ 


n 

^ a^Je A 


-1 


B 

(B2 + 4a2)1/2 




(A8) 


(A9) 

(AlO) 
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where 

_ PTC 
D = e' b 


a = 


Ku, 


b y'Tv/p^.^ 




u(l - u) 


/o (1 + Bu - A^O^) 3/2 

For a flat plate 


^s±D~^ B 

eA (b2 + 4a 2)1/2 du 


Tw Cf _ de 
^ dx 


Poo^« 


Therefore, from eguation (AlO) 


21= A. 

2 dx 


t*" R ^ 


B 


KP b 


or, from equation (A12), 


(b2+ 4a2)1/2 


B 


e = -ill- ga^JeA®^ (b2 + 4a2)^/2 

Kp„u« b 


(All) 


(A12) 


(A13) 


(A14) 


V/hen Rq is defined as 


Rfl = 


"oaP»e 


(AI 5 ) 


and 


equation (Al4) becomes 


Po. 


(g)- 




(A 16 ) 


Van Driest has shown that a first approximation to J obtained by 
integration by parts results in 


a2j = 


— sin 


- 1 2A^ - B 


(1 + B - A2)1/2 


(b2+ 4a2) 1/2 


(AIT) 
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Thus K 0 can be wrl'tten 



1 exp [ ^(sin"^cp 

(1 + B - V2 


+ sin ^ ijr) ] 


(AI 8 ) 


where 

2A^ - B ^ B 

(b^ + (B^ + 


Por the case of M«= 0 and Tw/Too = 1, equation (AI 8 ) was inte- 
grated to yield a relation between the average skin-friction coefficient 
and the Reynolds number based on the distance along the plate. Agree- 
ment with low-speed experl mental results is obtained I'Then F = 6.5* 

To determine how much the deviation of Pr and a from unity 
influences the local skin-friction coefficient, the results expressed 
by equation (AI 8 ) are shoivn in figures 3 and 4. The dashed lines repre- 
sent the relationship between cf/2 and Rq when Rr = 0»72 and r = O .89 
(a computed to agree with recovery factor). The solid lines represent 
the corresponding relationship for Rr = 1, a = 1 ana, conseq.uently, 
r - 1. It is observed that the deviation of Rr and a from unity pro- 
duces very little effect on Cf/2 at the lower Ifech numbers. At the 
higher Mach numbers the effect has some Reynolds number dependence; 
however, the largest deviation, between the two results is of the order 
of 7 percent. It should be noted that the ¥an Driest theory is kno>m 
to yield values of qkin-frictlon coefficient which are hi^er than 
measured values on the order of 11 percent at M = 2.5 (ref. 15)* The 
allowed variation of Rr or a does not iroprove this situation. 
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Figure /. — Modified Reynolds analogy for compressible boundary 
layer (insulated surface temperature). 
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Figure 2. — Modified Reynolds analogy for compressible boundary 
layer (wall temperature equal to free-stream temperature). 










